Comparison with Kapsenberg et al. (2017)

Description

Time series trends analysis by residuals (substracting climatological monthly means)

Kapsenberg et al. (2017)

2007-2015 : Comparison with Kapsenberg results

SURFACE

Table of monthly means (2007-2015) - surface data

➝ Comparaison supplement, Tab S1 (Kapsenberg) same period : results ok ✓

monthly_means_0m_0715 <- ungroup(data_0715) %>% 
  group_by(monthd) %>%
  summarise(
    Salinity_month = mean(salinity, na.rm = TRUE),
    Temperature_month = mean(temperature, na.rm = TRUE),
    dic_month = mean(dic, na.rm = TRUE),
    ta_month = mean(ta, na.rm = TRUE),
    pH_calc_month = mean(pH_calc, na.rm = TRUE),
    pH18_month = mean(pH18, na.rm = TRUE),
    pCO2_month = mean(pCO2, na.rm = TRUE),
    OmegaCalcite_month = mean(OmegaCalcite, na.rm = TRUE),
    OmegaAragonite_month = mean(OmegaAragonite, na.rm = TRUE),
    Nta_month = mean(Nta, na.rm = TRUE),
    Ndic_month = mean(dic, na.rm = TRUE))
monthly_means_0m_0715
## # A tibble: 12 × 12
##    monthd Salinity_month Temperature_month dic_month ta_month pH_calc_month
##     <dbl>          <dbl>             <dbl>     <dbl>    <dbl>         <dbl>
##  1      1           37.9              14.4     2251.    2551.          8.13
##  2      2           37.9              13.6     2260.    2556.          8.14
##  3      3           37.9              13.7     2260.    2555.          8.13
##  4      4           37.8              15.0     2254.    2553.          8.12
##  5      5           37.6              18.1     2241.    2548.          8.09
##  6      6           37.7              21.3     2233.    2544.          8.05
##  7      7           37.9              24.1     2228.    2549.          8.02
##  8      8           38.1              24.8     2227.    2556.          8.02
##  9      9           38.2              23.5     2230.    2562.          8.04
## 10     10           38.2              20.8     2227.    2561.          8.08
## 11     11           38.0              18.2     2232.    2555.          8.11
## 12     12           38.0              16.1     2242.    2556.          8.12
## # ℹ 6 more variables: pH18_month <dbl>, pCO2_month <dbl>,
## #   OmegaCalcite_month <dbl>, OmegaAragonite_month <dbl>, Nta_month <dbl>,
## #   Ndic_month <dbl>

Table of time series anomalies regression analyses (2007-2015) - surface data

➝ Comparison Table 2 (Kapsenberg) same period : results ok ✓

ano_0m_0715 <- left_join(ungroup(data_0715), monthly_means_0m_0715, by = "monthd") %>%
  mutate(
    Salinity_ano = salinity - Salinity_month,
    Temperature_ano = temperature - Temperature_month,
    dic_ano = dic - dic_month,
    ta_ano = ta - ta_month,
    pH_calc_ano = pH_calc - pH_calc_month,
    pH18_ano = pH18 - pH18_month,
    pCO2_ano = pCO2 - pCO2_month,
    OmegaCalcite_ano = OmegaCalcite - OmegaCalcite_month,
    OmegaAragonite_ano = OmegaAragonite - OmegaAragonite_month,
    Nta_ano = Nta - Nta_month,
    Ndic_ano = dic - Ndic_month)

# Regressions and tables of key parameters
var_list <-
  c("salinity",
    "temperature",
    "dic",
    "ta",
    "pH_calc",
    "pH18",
    "pCO2",
    "OmegaCalcite",
    "OmegaAragonite",
    "Nta",
    "Ndic")

lm.test <- vector("list", length(var_list))
 
 for(i in seq_along(var_list)){
     lm.test[[i]] <- lm(reformulate(var_list[i], "sampling_date"), data = ano_0m_0715)
 }


lms <- lapply(var_list, function(x) {
  summary(lm(substitute(i ~ decimal_date(sampling_date), list(i = as.name(x))), 
             data = ano_0m_0715))
})

reg <- NULL
for (i in 1:length(var_list)) {
  #one loops through all anomalies
  # calculate probability of fstatistic because it cannot be extracted from lms above
  # see http://stats.stackexchange.com/questions/92824/extracting-the-model-p-value-for-a-multiple-regression-in-r
  # slope returned by lm is per second
  prob <-
    pf(lms[[i]]$fstatistic[1],
       lms[[i]]$fstatistic[2],
       lms[[i]]$fstatistic[3],
       lower.tail = FALSE)
  reg <-
    rbind(reg, as.numeric(
      c(
        lms[[i]]$coefficients[2, 1],
        lms[[i]]$coefficients[2, 2],
        lms[[i]]$coefficients[2, 4],
        lms[[i]]$coefficients[1, 1],
        lms[[i]]$coefficients[1, 2],
        lms[[i]]$coefficients[1, 4],
        lms[[i]]$fstatistic[1],
        lms[[i]]$fstatistic[3],
        lms[[i]]$r.squared,
        prob
      )
    ))
}

colnames(reg) <-
  c("Slope",
    "SE Slope",
    "P Slope",
    "Intercept",
    "SE int.",
    "P int.",
    "F",
    "df",
    "R2",
    "P value")
row.names(reg) <- var_list
reg_ano_surf <- reg
slope_temp <- reg_ano_surf[2, 1]

# Regression and table anomalies
var_list <-
  c(
    "Salinity_ano",
    "Temperature_ano",
    "dic_ano",
    "ta_ano",
    "pH_calc_ano",
    "pH18_ano",
    "pCO2_ano",
    "OmegaCalcite_ano",
    "OmegaAragonite_ano",
    "Nta_ano",
    "Ndic_ano"
  )
lms <- lapply(var_list, function(x) {
  summary(lm(substitute(i ~ decimal_date(sampling_date), list(i = as.name(x))), 
             data = ano_0m_0715))
})
reg <- NULL
for (i in 1:length(var_list)) {
  #one loops through all anomalies
  # calculate probability of fstatistic because it cannot be extracted from lms above
  # see http://stats.stackexchange.com/questions/92824/extracting-the-model-p-value-for-a-multiple-regression-in-r
  # slope returned by lm is per second
  prob <-
    pf(lms[[i]]$fstatistic[1],
       lms[[i]]$fstatistic[2],
       lms[[i]]$fstatistic[3],
       lower.tail = FALSE)
  reg <-
    rbind(reg, as.numeric(
      c(
        lms[[i]]$coefficients[2, 1],
        lms[[i]]$coefficients[2, 2],
        lms[[i]]$coefficients[2, 4],
        lms[[i]]$coefficients[1, 1],
        lms[[i]]$coefficients[1, 2],
        lms[[i]]$coefficients[1, 4],
        lms[[i]]$fstatistic[1],
        lms[[i]]$fstatistic[3],
        lms[[i]]$r.squared,
        prob
      )
    ))
}
colnames(reg) <-
  c("Slope",
    "SE Slope",
    "P Slope",
    "Intercept",
    "SE int.",
    "P int.",
    "F",
    "df",
    "R2",
    "P value")
row.names(reg) <- var_list
reg_ano_surf <- reg
slope_temp <- reg_ano_surf[2, 1]
slope_pH <- reg_ano_surf[7, 1]

reg_ano_surf
##                           Slope     SE Slope      P Slope    Intercept
## Salinity_ano       -0.001689114 0.0043739500 6.995639e-01     3.397891
## Temperature_ano     0.071869400 0.0216700091 9.924272e-04  -144.576487
## dic_ano             3.159510272 0.1933457350 1.541892e-46 -6355.699197
## ta_ano              2.080353058 0.1867599667 2.197247e-25 -4184.893358
## pH_calc_ano        -0.002908521 0.0003337822 7.606349e-17     5.850844
## pH18_ano           -0.001884556 0.0002215758 3.540344e-16     3.791014
## pCO2_ano            3.627942537 0.3896601695 8.181341e-19 -7298.048505
## OmegaCalcite_ano   -0.012171909 0.0022940518 1.850430e-07    24.485279
## OmegaAragonite_ano -0.007271327 0.0015854405 6.022437e-06    14.627159
## Nta_ano             2.213010149 0.2738614887 7.131981e-15 -4451.749878
## Ndic_ano            3.159510272 0.1933457350 1.541892e-46 -6355.699197
##                        SE int.       P int.           F  df           R2
## Salinity_ano         8.7988267 6.995641e-01   0.1491318 416 0.0003583613
## Temperature_ano     43.5926333 9.924367e-04  10.9994216 411 0.0260650158
## dic_ano            388.9363365 1.542099e-46 267.0363299 408 0.3955880863
## ta_ano             375.6916128 2.197423e-25 124.0813516 414 0.2305996133
## pH_calc_ano          0.6714441 7.606755e-17  75.9308053 405 0.1578830146
## pH18_ano             0.4457270 3.540525e-16  72.3391874 405 0.1515467185
## pCO2_ano           783.8495602 8.181828e-19  86.6859612 405 0.1763035108
## OmegaCalcite_ano     4.6147686 1.850472e-07  28.1521109 405 0.0649935904
## OmegaAragonite_ano   3.1893094 6.022540e-06  21.0342677 405 0.0493722437
## Nta_ano            550.9074895 7.132321e-15  65.2987521 414 0.1362381017
## Ndic_ano           388.9363365 1.542099e-46 267.0363299 408 0.3955880863
##                         P value
## Salinity_ano       6.995639e-01
## Temperature_ano    9.924272e-04
## dic_ano            1.541892e-46
## ta_ano             2.197247e-25
## pH_calc_ano        7.606349e-17
## pH18_ano           3.540344e-16
## pCO2_ano           8.181341e-19
## OmegaCalcite_ano   1.850430e-07
## OmegaAragonite_ano 6.022437e-06
## Nta_ano            7.131981e-15
## Ndic_ano           1.541892e-46

Plot of anomalies : temperature (2007-2015) - surface

Plot of anomalies : salinity (2007-2015) - surface

Plot of anomalies : pH (2007-2015) - surface

Plot of anomalies : alkalinity (2007-2015) - surface

Plot of anomalies : pCO2 (2007-2015) - surface

DEPTH 50 m

Table of monthly means (2007-2015) - 50m data

➝ Comparaison supplement, Tab S1 (Kapsenberg) same period : results ok ✓

## # A tibble: 12 × 12
##    monthd Salinity_month Temperature_month dic_month ta_month pH_calc_month
##     <dbl>          <dbl>             <dbl>     <dbl>    <dbl>         <dbl>
##  1      1           38.0              14.5     2253.    2555.          8.13
##  2      2           38.0              13.7     2260.    2556.          8.13
##  3      3           38.0              13.6     2261.    2555.          8.13
##  4      4           38.0              14.1     2258.    2553.          8.12
##  5      5           38.0              15.0     2253.    2551.          8.12
##  6      6           38.0              15.1     2249.    2549.          8.12
##  7      7           38.0              15.4     2246.    2546.          8.11
##  8      8           38.0              15.3     2241.    2544.          8.12
##  9      9           38.0              16.2     2237.    2548.          8.12
## 10     10           38.1              18.5     2233.    2555.          8.10
## 11     11           38.1              18.0     2234.    2557.          8.11
## 12     12           38.0              16.2     2242.    2557.          8.12
## # ℹ 6 more variables: pH18_month <dbl>, pCO2_month <dbl>,
## #   OmegaCalcite_month <dbl>, OmegaAragonite_month <dbl>, Nta_month <dbl>,
## #   Ndic_month <dbl>

Table of time series anomalies regression analyses (2007-2015) - 50m data

➝ Comparison Table 2 (Kapsenberg) same period: results ok ✓

##                           Slope     SE Slope      P Slope    Intercept
## Salinity_ano        0.006280090 0.0020205419 2.013976e-03   -12.633358
## Temperature_ano     0.087907062 0.0187730250 3.867940e-06  -176.840381
## dic_ano             2.237743420 0.2042325862 1.192951e-24 -4501.572889
## ta_ano              1.612015807 0.1533532351 4.883020e-23 -3242.803820
## pH_calc_ano        -0.002647707 0.0002445368 3.749752e-24     5.326314
## pH18_ano           -0.001305946 0.0002666917 1.412206e-06     2.627133
## pCO2_ano            2.813812116 0.2453557568 1.528315e-26 -5660.463142
## OmegaCalcite_ano   -0.007213340 0.0026443975 6.654962e-03    14.510865
## OmegaAragonite_ano -0.003929317 0.0018128827 3.078619e-02     7.904491
## Nta_ano             1.187830184 0.1244390877 1.254187e-19 -2389.492859
## Ndic_ano            2.237743420 0.2042325862 1.192951e-24 -4501.572889
##                        SE int.       P int.          F  df         R2
## Salinity_ano         4.0646310 2.013993e-03   9.660421 410 0.02301962
## Temperature_ano     37.7652440 3.868010e-06  21.926990 406 0.05124003
## dic_ano            410.8463170 1.193045e-24 120.052305 406 0.22821363
## ta_ano             308.4925535 4.883378e-23 110.497565 409 0.21270083
## pH_calc_ano          0.4919279 3.750038e-24 117.233404 403 0.22534771
## pH18_ano             0.5364963 1.412234e-06  23.979014 403 0.05615970
## pCO2_ano           493.5753768 1.528442e-26 131.521782 403 0.24605505
## OmegaCalcite_ano     5.3196612 6.655007e-03   7.440796 403 0.01812879
## OmegaAragonite_ano   3.6469260 3.078633e-02   4.697802 403 0.01152275
## Nta_ano            250.3274998 1.254266e-19  91.116090 409 0.18218988
## Ndic_ano           410.8463170 1.193045e-24 120.052305 406 0.22821363
##                         P value
## Salinity_ano       2.013976e-03
## Temperature_ano    3.867940e-06
## dic_ano            1.192951e-24
## ta_ano             4.883020e-23
## pH_calc_ano        3.749752e-24
## pH18_ano           1.412206e-06
## pCO2_ano           1.528315e-26
## OmegaCalcite_ano   6.654962e-03
## OmegaAragonite_ano 3.078619e-02
## Nta_ano            1.254187e-19
## Ndic_ano           1.192951e-24

Plot of anomalies : temperature (2007-2015) - 50m

Plot of anomalies : salinity (2007-2015) - 50m

Plot of anomalies : pH (2007-2015) - 50m

Plot of anomalies : alcalinity (2007-2015) - 50m

Plot of anomalies : pCO2 (2007-2015) - 50m

2007-2022

→ SOMLIT clean + interpolations ok ☑

SURFACE

Table of monthly means (2007-2022) - surface data

Table of time series anomaly regression analyses (2007-2022) - surface data

##                         Slope     SE Slope       P Slope     Intercept
## Temperature_ano  5.986065e-02 0.0090266801  5.975146e-11   -120.617872
## Salinity_ano     5.640707e-05 0.0020994838  9.785722e-01     -0.113659
## pH_calc_ano     -3.964373e-03 0.0001547736 9.561431e-107      7.987629
## pH18_ano        -3.113001e-03 0.0001180322 1.869102e-111      6.272240
## ta_ano           5.533254e-01 0.0961513994  1.226052e-08  -1114.869198
## dic_ano          2.439697e+00 0.1059150698  6.601288e-91  -4915.630099
## pCO2_ano         5.081188e+00 0.2104165592  6.711138e-95 -10238.096229
##                     SE int.        P int.            F  df           R2
## Temperature_ano  18.1886073  5.975814e-11 4.397711e+01 832 5.020349e-02
## Salinity_ano      4.2304243  9.785722e-01 7.218418e-04 832 8.675975e-07
## pH_calc_ano       0.3118469 9.570277e-107 6.560775e+02 819 4.447749e-01
## pH18_ano          0.2378181 1.870888e-111 6.955968e+02 819 4.592620e-01
## ta_ano          193.7314049  1.226154e-08 3.311692e+01 819 3.886429e-02
## dic_ano         213.4038134  6.606689e-91 5.305866e+02 819 3.931475e-01
## pCO2_ano        423.9880107  6.810126e-95 5.831370e+02 722 4.468014e-01
##                       P value
## Temperature_ano  5.975146e-11
## Salinity_ano     9.785722e-01
## pH_calc_ano     9.561431e-107
## pH18_ano        1.869102e-111
## ta_ano           1.226052e-08
## dic_ano          6.601288e-91
## pCO2_ano         6.711138e-95

Plot of anomalies : Temperature (2007-2022) - surface

Plot of anomalies : salinity (2007-2022) - surface

Plot of anomalies : pH_calc (2007-2022) - surface

Plot of anomalies : pH(18) (2007-2022) - surface

Plot of anomalies : alkalinity (2007-2022) - surface

Plot of anomalies : DIC (2007-2022) - surface

Plot of anomalies : pCO2 (2007-2022) - surface

DEPTH 50 m

Table of monthly means (2007-2022) - 50m data

(rajouter les SD)

## # A tibble: 12 × 15
##    Month Salinity_month SD_sal Temperature_month SD_temp dic_month SD_dic
##    <dbl>          <dbl>  <dbl>             <dbl>   <dbl>     <dbl>  <dbl>
##  1     1           38.1 0.123               14.6   0.863     2259.   18.4
##  2     2           38.1 0.121               13.8   0.388     2271.   17.1
##  3     3           38.1 0.134               13.7   0.387     2271.   15.7
##  4     4           38.0 0.113               14.2   0.521     2268.   17.6
##  5     5           38.0 0.120               15.1   0.806     2263.   19.7
##  6     6           38.0 0.0975              15.2   0.672     2258.   18.1
##  7     7           38.0 0.0826              15.4   0.920     2256.   19.0
##  8     8           38.0 0.0743              15.4   0.590     2253.   16.9
##  9     9           38.0 0.103               16.2   1.71      2246.   18.4
## 10    10           38.1 0.170               18.5   1.76      2236.   14.9
## 11    11           38.1 0.124               18.1   1.13      2237.   14.9
## 12    12           38.1 0.140               16.1   0.968     2249.   14.2
## # ℹ 8 more variables: ta_month <dbl>, SD_ta <dbl>, pH_calc_month <dbl>,
## #   SD_PH <dbl>, pH18_month <dbl>, SD_PH18 <dbl>, pCO2_month <dbl>,
## #   SD_pCO2 <dbl>

Table of time series anomaly regression analyses (2007-2022) - 50m data

##                        Slope     SE Slope       P Slope    Intercept
## Temperature_ano  0.026456564 0.0073678259  3.489534e-04   -53.309392
## Salinity_ano     0.009666510 0.0008266589  2.362696e-29   -19.477804
## pH_calc_ano     -0.003343820 0.0001278701 4.579139e-110     6.737306
## pH18_ano        -0.002980286 0.0001380500  3.547310e-82     6.004839
## ta_ano           0.583934397 0.0803203798  8.396021e-13 -1176.541935
## dic_ano          2.289966337 0.1044826549  3.779468e-84 -4613.945403
## pCO2_ano         3.454482642 0.1256618844 2.184751e-118 -6960.274502
##                     SE int.        P int.         F  df         R2
## Temperature_ano  14.8460452  3.489658e-04  12.89401 832 0.01526110
## Salinity_ano      1.6657037  2.363428e-29 136.73722 832 0.14114996
## pH_calc_ano       0.2576402 4.583474e-110 683.82918 819 0.45502788
## pH18_ano          0.2781511  3.549988e-82 466.06148 819 0.36267641
## ta_ano          161.8341636  8.397098e-13  52.85385 819 0.06062237
## dic_ano         210.5177179  3.782376e-84 480.36328 819 0.36969128
## pCO2_ano        253.1908589 2.186932e-118 755.71648 819 0.47990638
##                       P value
## Temperature_ano  3.489534e-04
## Salinity_ano     2.362696e-29
## pH_calc_ano     4.579139e-110
## pH18_ano         3.547310e-82
## ta_ano           8.396021e-13
## dic_ano          3.779468e-84
## pCO2_ano        2.184751e-118

Plot of anomalies : Temperature (2007-2022) - 50m

Plot of anomalies : salinity (2007-2022) - 50m

Plot of anomalies : pH_calc (2007-2022) - 50m

Plot of anomalies : pH(18) (2007-2022) - 50m

Plot of anomalies : alkalinity (2007-2022) - 50m

Plot of anomalies : DIC (2007-2022) - 50m

Plot of anomalies : pCO2 (2007-2022) - 50m

1995-2022 : temperature & salinity

05-92 to 12-22 but not enough data before Jan. 1995

SURFACE

Table of monthly means for salinity and temperature (1995-2022) - 0m data

## # A tibble: 12 × 5
##    month Salinity_month SD_sal Temperature_month SD_temp
##    <chr>          <dbl>  <dbl>             <dbl>   <dbl>
##  1 01              38.0  0.288              14.3   0.699
##  2 02              38.0  0.192              13.6   0.419
##  3 03              38.0  0.201              13.6   0.498
##  4 04              37.9  0.238              14.8   0.934
##  5 05              37.7  0.255              17.8   1.54 
##  6 06              37.7  0.277              21.6   1.88 
##  7 07              37.9  0.302              24.4   1.94 
##  8 08              38.1  0.331              25.3   1.33 
##  9 09              38.2  0.176              23.4   1.39 
## 10 10              38.1  0.214              20.8   1.16 
## 11 11              38.0  0.370              18.1   1.07 
## 12 12              38.0  0.329              16.0   0.951

Table of salinity and temperature time series anomaly regression analyses (1995-2022) - 0m data

##                        Slope     SE Slope      P Slope  Intercept  SE int.
## Salinity_ano    -0.001033649 0.0009029316 2.525074e-01   2.076398 1.813827
## Temperature_ano  0.038182803 0.0040684385 2.570245e-20 -76.701761 8.172761
##                       P int.         F   df           R2      P value
## Salinity_ano    2.525113e-01  1.310498 1350 0.0009697979 2.525074e-01
## Temperature_ano 2.571986e-20 88.080569 1350 0.0612487026 2.570245e-20

Plot of anomalies : Temperature (1995-2022) - 0m

Plot of anomalies : salinity (1995-2022) - 0m

DEPTH 50

Table of monthly means for salinity and temperature (1995-2022) - 50m data

## # A tibble: 12 × 5
##    month Salinity_month SD_sal Temperature_month SD_temp
##    <chr>          <dbl>  <dbl>             <dbl>   <dbl>
##  1 01              38.1 0.130               14.3   0.668
##  2 02              38.1 0.121               13.6   0.417
##  3 03              38.1 0.125               13.5   0.421
##  4 04              38.0 0.107               14.0   0.539
##  5 05              38.0 0.117               14.8   0.765
##  6 06              38.0 0.104               15.1   0.635
##  7 07              38.0 0.0982              15.2   0.633
##  8 08              38.0 0.101               15.4   0.689
##  9 09              38.0 0.110               16.1   1.54 
## 10 10              38.1 0.165               18.6   1.81 
## 11 11              38.1 0.129               17.9   1.14 
## 12 12              38.1 0.132               16.0   0.953

Table of salinity and temperature time series anomaly regression analyses (1995-2022) - 50m data

##                       Slope     SE Slope      P Slope  Intercept   SE int.
## Salinity_ano    0.003886329 0.0003942302 3.490107e-22  -7.806559 0.7919055
## Temperature_ano 0.019299605 0.0031399001 1.040676e-09 -38.767620 6.3072498
##                      P int.        F   df         R2      P value
## Salinity_ano    3.49267e-22 97.18058 1347 0.06729115 3.490107e-22
## Temperature_ano 1.04099e-09 37.78028 1349 0.02724316 1.040676e-09

Plot of anomalies : Temperature (1995-2022) - 50m

Plot of anomalies : salinity (1995-2022) - 50m

2000-2022 : Analysis every 10 years

SURFACE

Jan.2000 - Dec.2009

Table of monthly means (2000-2009) - 0m data

## # A tibble: 12 × 5
##    month Salinity_month SD_sal Temperature_month SD_temp
##    <chr>          <dbl>  <dbl>             <dbl>   <dbl>
##  1 01              38.1  0.123              14.1   0.643
##  2 02              38.1  0.131              13.4   0.460
##  3 03              38.1  0.107              13.6   0.525
##  4 04              37.9  0.156              14.5   0.939
##  5 05              37.7  0.290              18.0   1.50 
##  6 06              37.8  0.284              21.7   1.87 
##  7 07              38.0  0.217              23.9   1.83 
##  8 08              38.2  0.193              25.0   1.34 
##  9 09              38.2  0.161              23.1   1.15 
## 10 10              38.1  0.209              20.8   1.16 
## 11 11              38.0  0.200              18.1   1.03 
## 12 12              38.0  0.178              16.1   0.845

Table of anomaly regression analyses (2000-2009) - 0m data

##                       Slope    SE Slope    P Slope  Intercept   SE int.
## Salinity_ano    -0.00188507 0.003123799 0.54649521   3.779479  6.263081
## Temperature_ano  0.03810298 0.018994081 0.04541915 -76.394730 38.082303
##                     P int.        F  df           R2    P value
## Salinity_ano    0.54649562 0.364157 474 0.0007676739 0.54649521
## Temperature_ano 0.04541937 4.024216 474 0.0084184354 0.04541915

Plot of anomalies : Temperature (2000-2009) - 0m

Plot of anomalies : Salinity (2000-2009) - 0m

Jan.2010 - Dec.2022

Table of monthly means (2010-2022) - 0m data

## # A tibble: 12 × 5
##    month Salinity_month SD_sal Temperature_month SD_temp
##    <chr>          <dbl>  <dbl>             <dbl>   <dbl>
##  1 01              38.0  0.375              14.5   0.711
##  2 02              38.0  0.233              13.7   0.375
##  3 03              37.9  0.245              13.7   0.482
##  4 04              37.8  0.295              15.1   0.927
##  5 05              37.8  0.254              17.7   1.59 
##  6 06              37.7  0.296              21.8   2.02 
##  7 07              37.8  0.373              25.0   1.76 
##  8 08              38.0  0.450              25.6   1.45 
##  9 09              38.2  0.180              23.7   1.53 
## 10 10              38.2  0.179              20.9   1.10 
## 11 11              37.9  0.518              18.3   1.00 
## 12 12              38.0  0.458              16.3   0.949

Table of anomaly regression analyses (2010-2022) - 0m data

##                       Slope    SE Slope      P Slope   Intercept   SE int.
## Salinity_ano    0.003813676 0.003557838 2.841797e-01   -7.689935  7.174074
## Temperature_ano 0.092233433 0.013161962 6.362107e-12 -185.980424 26.539958
##                       P int.         F  df          R2      P value
## Salinity_ano    2.841805e-01  1.148987 618 0.001855753 2.841797e-01
## Temperature_ano 6.362608e-12 49.106106 618 0.073610638 6.362107e-12

Plot of anomalies : Temperature (2010-2022) - 0m

Plot of anomalies : Salinity (2010-2022) - 0m

DEPTH 50

Jan.2000 - Dec.2009

Table of monthly means (2000-2009) - 50m data

## # A tibble: 12 × 5
##    month Salinity_month SD_sal Temperature_month SD_temp
##    <chr>          <dbl>  <dbl>             <dbl>   <dbl>
##  1 01              38.1 0.0914              14.2   0.654
##  2 02              38.1 0.0779              13.5   0.461
##  3 03              38.1 0.0799              13.4   0.473
##  4 04              38.1 0.0846              13.8   0.528
##  5 05              38.0 0.118               14.8   0.859
##  6 06              38.0 0.111               14.9   0.622
##  7 07              38.0 0.113               14.9   0.613
##  8 08              38.0 0.137               15.2   0.528
##  9 09              38.0 0.102               15.8   1.10 
## 10 10              38.1 0.176               18.3   1.92 
## 11 11              38.1 0.144               18.0   1.02 
## 12 12              38.0 0.131               16.0   0.813

Table of anomaly regression analyses (2000-2009) - 50m data

##                        Slope    SE Slope    P Slope  Intercept   SE int.
## Salinity_ano    -0.002918298 0.001847716 0.11490525   5.851096  3.704616
## Temperature_ano  0.027183103 0.013892758 0.05097301 -54.501262 27.854567
##                     P int.        F  df          R2    P value
## Salinity_ano    0.11490562 2.494531 477 0.005202419 0.11490525
## Temperature_ano 0.05097325 3.828434 477 0.007962162 0.05097301

Plot of anomalies : Temperature (2000-2009) - 50m

Plot of anomalies : Salinity (2000-2009) - 50m

Jan.2010 - Dec.2022

Table of monthly means (2010-2022) - 50m data

## # A tibble: 12 × 5
##    month Salinity_month  SD_sal Temperature_month SD_temp
##    <chr>          <dbl>   <dbl>             <dbl>   <dbl>
##  1 01              38.1  0.134               14.5   0.663
##  2 02              38.1  0.128               13.8   0.356
##  3 03              38.1  0.146               13.6   0.349
##  4 04              38.0  0.124               14.2   0.517
##  5 05              38.1  0.116               15.0   0.703
##  6 06              38.1  0.0869              15.1   0.607
##  7 07              38.1  0.0752              15.3   0.643
##  8 08              38.0  0.0593              15.4   0.640
##  9 09              38.1  0.0984              16.2   1.92 
## 10 10              38.1  0.152               18.9   1.81 
## 11 11              NA   NA                   18.0   1.19 
## 12 12              NA   NA                   16.3   1.00

Table of anomaly regression analyses (2010-2022) - 50m data

##                     Slope    SE Slope      P Slope Intercept   SE int.
## Salinity_ano    0.0114991 0.001276746 3.963691e-18 -23.18538  2.574279
## Temperature_ano 0.0257115 0.010697870 1.653841e-02 -51.84285 21.570465
##                       P int.         F  df          R2      P value
## Salinity_ano    3.964160e-18 81.118272 523 0.134275482 3.963691e-18
## Temperature_ano 1.653859e-02  5.776438 613 0.009335259 1.653841e-02

Plot of anomalies : Temperature (2010-2022) - 50m

Plot of anomalies : Salinity (2010-2022) - 50m


Comparison SOMLIT - EOL results : 2018-2022

SURFACE

temperature (°C) & salinity

Table and plot of monthly means (2017-2022) - EOL temperature & salinity

monthly_means_EOL_1822 <- ungroup(EOL_hourly_1822) %>% 
  group_by(Month) %>%
  summarise(
    Salinity_month = mean(salinity, na.rm = TRUE),
    SD_sal = sd(salinity, na.rm = TRUE),
    Temperature_month = mean(temperature, na.rm = TRUE),
    SD_temp = sd(temperature, na.rm = TRUE))

monthly_means_EOL_1822
## # A tibble: 12 × 5
##    Month Salinity_month SD_sal Temperature_month SD_temp
##    <chr>          <dbl>  <dbl>             <dbl>   <dbl>
##  1 01              38.0  0.107              14.6   0.564
##  2 02              38.1  0.151              13.9   0.277
##  3 03              38.0  0.268              13.9   0.388
##  4 04              37.9  0.303              15.1   1.08 
##  5 05              37.8  0.202              18.2   1.92 
##  6 06              37.8  0.231              22.4   1.80 
##  7 07              37.9  0.150              26.1   1.40 
##  8 08              38.1  0.195              26.7   1.01 
##  9 09              38.2  0.167              24.5   1.44 
## 10 10              38.1  0.182              21.2   1.10 
## 11 11              37.7  3.13               18.7   1.05 
## 12 12              38.0  0.160              16.4   0.765
monthly_means_somlit_1822 <- ungroup(SOMLIT_1822) %>%
  group_by(Month) %>%
  summarise(
    Salinity_month = mean(impute_sal, na.rm = TRUE),
    SD_sal = sd(impute_sal, na.rm = TRUE),
    Temperature_month = mean(impute_temp, na.rm = TRUE),
    SD_temp = sd(impute_temp, na.rm = TRUE))

monthly_means_somlit_1822
## # A tibble: 12 × 5
##    Month Salinity_month SD_sal Temperature_month SD_temp
##    <chr>          <dbl>  <dbl>             <dbl>   <dbl>
##  1 01              38.1  0.118              14.6   0.600
##  2 02              38.1  0.124              13.8   0.310
##  3 03              38.0  0.194              13.8   0.405
##  4 04              38.0  0.151              14.8   0.718
##  5 05              37.7  0.296              17.9   1.91 
##  6 06              37.5  0.548              22.4   1.88 
##  7 07              37.6  0.426              26.1   1.25 
##  8 08              37.8  0.556              26.2   1.51 
##  9 09              38.1  0.220              24.4   1.35 
## 10 10              38.1  0.238              21.2   1.06 
## 11 11              38.0  0.418              18.9   1.02 
## 12 12              38.0  0.191              16.3   0.754

Table of anomalies regression analyses (2018-2022) - EOL temperature & salinity

##                       Slope    SE Slope      P Slope  Intercept  SE int.
## Salinity_ano    0.002164672 0.003343598 5.173718e-01  -4.373596 6.755563
## Temperature_ano 0.031911638 0.004104452 7.733525e-15 -64.477267 8.293022
##                       P int.         F    df           R2      P value
## Salinity_ano    5.173726e-01  0.419137 38850 1.078848e-05 5.173718e-01
## Temperature_ano 7.733608e-15 60.448817 39611 1.523736e-03 7.733525e-15

Table of anomalies regression analyses (2018-2022) - SOMLIT temperature & salinity

##                        Slope   SE Slope   P Slope  Intercept   SE int.
## Salinity_ano    -0.004897158 0.01386080 0.7241470   9.894662  28.00563
## Temperature_ano  0.039645714 0.04992486 0.4278683 -80.103804 100.87275
##                    P int.         F  df          R2   P value
## Salinity_ano    0.7241470 0.1248278 256 0.000487371 0.7241470
## Temperature_ano 0.4278684 0.6306070 256 0.002457255 0.4278683

EOL : Plot of anomalies : Temperature (°C) (2018-2022)

EOL : Plot of anomalies : Salinity (2018-2022)

without interpolation

SOMLIT : Plot of anomalies : Temperature (°C) (2018-2022)

SOMLIT : Plot of anomalies : Salinity (2018-2022)